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Abstract 

We investigate the effect of unparticle self- interactions at the Large Hadron Collider (LHC). 
Especially, we discuss the three-point correlation function, which is determined by conformal sym- 
metry up to a constant, and study its relation to processes with four-particle final states. These 
processes could be used as a favorable way to look for unparticle physics, and for weak enough 
couplings to the Standard Model, even the only way. We find updated upper bounds on the cross 
sections for unparticle-mediated 47 final states at the LHC and novel upper bounds for the corre- 
sponding 272^ and Ai final states. The size of the allowed cross sections obtained are comparably 
large for large values of the scaling dimension of the unparticle sector, but they decrease with 
decreasing values of this parameter. In addition, we present relevant distributions for the different 
final states, enabling the possible identification of the unparticle scaling dimension if there was to 
be a large number of events of such final states at the LHC. 
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I. INTRODUCTION 



The Standard Model (SM) of particle physics has to date been extremely successful 
yielding consistent results and most accurate predicting properties of particles and their 
interactions. However, in its minimal version, it does not contain massive neutrinos, gravity, 
and/or dark matter. Although it should be fairly easy to include massive neutrinos as 
well as neutrino oscillations in a minimally extended version of the SM, one of the most 
important particles of the SM has not yet been experimentally confirmed, namely the Higgs 
boson, which is responsible for giving masses to the particles of the SM. Therefore, the Large 
Hadron Collider (LHC) has been constructed and will hopefully be able to discover it. In the 
connection with the process of constructing the LHC, a lot of alternative pictures have been 
proposed, including for example unparticle physics that was introduced by Howard Georgi in 
2007 [1] . Thus, in the last couple of years, several studies of the phenomenology of unparticle 
physics have appeared in the literature. In this work, we will focus on studying the effects of 
unparticle self-interactions at colliders, and especially the LHC. In general, our investigation 
will be based on the unparticle three-point correlation function, and in particular, we will 
discuss the phenomenology of unparticle self-interactions in processes with different four- 
particle final states such as four photons (47), two photons and two charged Icptons (272^), 
and four charged leptons (4£) as well as final states including neutrinos, i.e., 2'-)vV and 2£z/z7. 

Previously, an investigation of unparticle self-interactions that is of great importance for 
this work has been presented in the following paper: In Ref. [2], collider signals of unparticle 
self-interactions especially with four-photon final states are investigated, which lead to upper 
bounds on cross sections using Tevatron data that follow from restricting the unparticle self- 
interactions. In addition, Strassler [3] has discussed self-interactions in general and Georgi 
and Katz [4] have discussed the production of unparticles through the amputated three-point 
function. 

Our work is organized as follows. In Sec. II, we introduce our model based on a conformal 
hidden sector coupled to the SM and also present the specific couphngs to the SM. Then, 
in Sec. Ill, we explicitly derive the unparticle three-point correlation function and discuss 
the general form of the four-point correlation function. Next, in Sec. IV, we present the 
phenomenology of our model concerning different four-particle final states mediated by the 
unparticle three-point function at the LHC. By using existing Tevatron data, we constrain 
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the constant appearing in the three-point function, enabhng determination of the maximum 
allowed cross sections at the LHC. Furthermore, relevant distributions are given, which could 
be used to identify unparticle physics as a possible interpretation of future data in the event 
that a large number of events with four-body final states were to be measured at the LHC. 
Finally, in Sec. V, we give a short summary of our work as well as our main conclusions. 



II. THE MODEL 



The most famous examples of four-dimensional conformal field theories, which could serve 
as hidden sectors, are theories of so-called Banks-Zaks type [5]. These are gauge theories 
with the number of fermions chosen such that an infrared stable fixed point emerges at two 
loops. In addition, one can have supersymmetric versions of these types of theories and even 
calculate the fraction of them being conformal in the infrared regime [6]. 

In general, unitary four-dimensional scale-invariant quantum field theories are also con- 
formally invariant with no known counterexamples [7-9] . Unitarity demands that the scaling 
dimension of an operator with Lorentz spins ji and j2 is d > ji -l- j2 + 2 — Sj^^j^^ [10]. This 
translates into ds > 1, dsp > 3/2, and > 3 for Lorentz scalar, spinor, and gauge invariant 
vector operators, respectively. However, much of the phenomenology of unparticle physics 
dealing with vector operators are performed for values violating these bounds. One also 
finds that the fields satisfy equations of motion for free fields and conservation laws if and 
only if the lower bounds are saturated. For vector operators, this includes the existence of 
zero-norm descendants = 0, which does not hold away from the unitarity bound. 

The basic scheme for our model is the following. There is a hidden sector of Banks-Zaks 
type that couples to the SM through a 'messenger field' with a large mass M, which simply 
means that both the SM and the hidden sector couple to this field. Below the scale M, one 
can use effective field theory and integrate out the heavy field, whereby one ends up with 
effective operators suppressed by powers of M on the form^ 

<j^n+duv-4^nOuV, (1) 

where are dimensionless constants and and Ouv are local operators constructed out of 
the SM and hidden sector fields with mass dimensions n and duv, respectively. The hidden 



^ Our notation will follow that of Ref. [Ill 
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sector becomes conformal at the scale Ku,^ for example, by quantum effects via dimensional 
transmutation. The operator Ouv changes into A^'^'''^Ou, where Ou has scaling and mass 
dimension d. This implies that the couplings in Eq. (1) have to be replaced by 

<j^n+duv-4^n(^^ = <^„+d_4^n^W' (2) 



n 



where the effective scale A„ is given by 

K^'-'=\^- (3) 

Note that M, Au, d, and duv are parameters of the messenger and hidden sector, while n 
will depend on what SM operator is under consideration. Thus, A„ will depend on n (and 
hence the subscript), i.e., operators of different mass dimension will couple with different 
strengths. In fact, these different energy scales counteract the effective field theory intuition 
that higher-dimensional operators are more highly suppressed [11]. In addition, note that 
the ehmination of the heavy field induces contact interactions between generic SM fields 
such as 

suppressed by powers of M. These interactions have the potential to drown any unparticle 
effects [7], at least for values of d much larger than 1, and thus could give the most stringent 
constraints on the parameter space relevant for unparticle physics. 

Since the scahng dimension of the scalar unparticle operator exceeds one and all the SM 
operators except one that it can couple to all have mass dimension n > 3, these couplings 
are all irrelevant. The only coupling that can possibly be relevant is the Higgs- unparticle 
coupling 

\duv-d -j 

For d < 2, this coupling is indeed relevant and can significantly alter the low-energy physics 
of the unparticle sector. Once the Higgs field develops a vacuum expectation value v (giving 
H'^H — > f^), this operator introduces a scale into the conformal field theory and will cause 
the theory to become nonconformal at the energy scale A^ [12], where 

, -l-d _ \ 2-d„.2 



Ay = (5) 



^ This scale can for example be defined as the scale, where the coupling has reached ^2/3 of its fixed-point 
value. 
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Below the scale A^^, the hidden sector presumably becomes a conventional particle sector. 
For consistency, one requires < Ku and the two scales to be reasonably well separated 
in order to give a window, where the sector is conformal. 

The standard way to treat the breaking of conformal symmetry is to introduce a mass 
gap in the spectral density. In the absence of fine tuning, this mass gap is at least several 
GeVs. This implies that there are no long-range forces from unparticle exchange and that 
precision constraints in practice disappear. Thus, low-energy experiments are not sensitive to 
unparticle effects, and these are most easily probed at coUiders. If the energy of the process 
in question is high enough, one can simply ignore the mass gap in the expression for the 
propagator. In Sec. IV, the unparticle three-point function will be used to study possible 
signals at the LHC Although it is unclear how to implement the breaking of conformal 
invariance into this function, it is a reasonable assumption that one can ignore the existence 
of the mass gap when calculating the three-point function as long as all the momentum 
invariants [g^ , gf|, and (?i+?2)^] in the function are much larger than the conformal symmetry 
breaking scale. Since this will generally be the case at the LHC, this subject will not be 
discussed any further. 

Covariance under conformal transformations determines the (exact) two-point functions 
of primary scalar operators Oi with scaling dimensions di to be 

mx)Om) « (Z^^ « / (^( V - iet-'e-^-^ (6) 

which give the expressions for the momentum space propagators that can enter in processes 
with internal unparticle lines. The form of the propagator can give unusual kinematical 
distributions due to the unusual dependence on and the unusual phase, i.e., e^''"^, which 
can lead to strange interference effects with photon and Z boson propagators [13-15]. Fur- 
thermore, the phase space is given by the imaginary part of the propagator 

The real proportionality constant (generally denoted by A^) corresponds to field normaliza- 
tion and is most commonly chosen such that the phase space interpolates that of d massless 
particles [1] . The unparticle phase space will enter in processes with an unparticle line in the 
final state. However, for d > 2, it is ultraviolet sensitive and one sometimes finds singular 
behavior. For example, the energy density at finite temperature and some cross sections are 
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proportional to 2 — d [16, 17]. This together with the lower bound on the scaling dimension 
from unitarity implies that one usually restricts the attention to the range 1 < d < 2 for 
scalars. 

Finally, it is generally assumed that an unparticle does not decay into SM particles. 
Thus, the unparticle will simply leave any detector undetected and lead to missing energy 
signals. However, the couplings to SM fields can through loops give a width to the unparticle 
propagator [18, 19], which means that the unparticle can, depending on model parameters, 
decay back into SM particles. For an overview of the phenomenology of unparticle physics, 
see for example Ref. [20] and references therein. 

Higher-order processes in the unparticle sector will enter through the unparticle n-point 
correlation function for n > 3. Some general ideas can be found in Ref. [3], and in Ref. [4], 
the production of unparticle stuff through its self-interactions has been studied. In this 
paper, we will follow Ref. [2] and look at the imphcations for LHC physics of the unparticle 
three-point function mediating interactions between SM particles. 



A. Couplings to the Standard Model 

In order to calculate obscrvablcs, one needs explicit expressions for the couplings of the 
unparticle operator to the SM. Initially, the only constraints on the SM operators arc that 
the interaction terms form Lorentz scalars. To reduce the number of possible interactions, 
one can assume that the unparticle operator is gauge invariant, and thus, the SM operators 
have to be so as well. All couplings of scalar, spinor, and vector unparticle operators to 
gauge invariant SM operators with dimensions less than or equal to four have been given in 
Ref. [21]. 

For a scalar unparticle field with scaling dimension d, the leading order (in A) interaction 
with SM gauge bosons is through [11] 

^F^^uF'^'Ou (8) 
and with fermions after electroweak symmetry breaking through 

-j^vfJnOu + -j^fL,RUL,Rd'Ou, (9) 

where the electromagnetic coupling has been pulled out for convenience. The first of these 
two operators is proportional to v. Integrating the second one by parts, the vector contri- 
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bution vanishes, while the axial- vector contribution is proportional to the fermion mass m/. 
Since nif <^ v ior all fermions except the top quark, the first interaction term will dominate 
for C4 = 0(c4^'^). Thus, the Feynman rules for the interactions of photons, gluons, and 
fermions with a scalar unparticle field are given by 

Ougg^Ou^^ vertex: . p,^/^- + p-p^), (10) 

/ 

OuJf vertex: ^^^Pr, (H) 

respectively, which will be used to study unparticle phenomenology in the rest of this work. 
For primary vectors, the interaction with SM fermions would be through 



C3 



fL,R 



lid- 

-'^3 



-rfL,RUL,RO^, (12) 



which would not be suppressed as the corresponding operator in Eq. (9). A vector unparticle 
field cannot couple directly to gauge bosons, since the number of uncontr acted Lorentz 
indices will always be odd. However, it can couple to them through a term hke H^D^HO^. A 
gauge invariant unparticle spinor can only couple in the same way as a right-handed neutrino, 
and thus, it is possible to replace the ordinary right-handed neutrino with a conformal one 
[22]. 



III. THE UNPARTICLE THREE-POINT CORRELATION FUNCTION 

In this section, we will investigate the scalar unparticle three-point correlation function. 
Then, in Sec. IV, we will make use of it to study unparticle self-interactions at the LHC. In 
general, the three-point function can involve fields with different spins [23] , but here we will 
only consider scalars. The three-point functions of scalar fields Oi with scaling dimensions 
di are constrained by conformal invariance as 

(C»i(Xi)C»2(x2)03(^3)) = (^2^)(<ii+d.-d3)/2(^2^)(<i.+d3-<ii)/2(^2^)(<ii+d3-<i.)/2 ' (l^) 

for some constant C123 and where Xij = Xi — xj. 

To calculate amplitudes, one needs the Fourier transform of the three-point function. 
This can be performed in a general number of dimensions D and for different values of the 
scaling dimensions di of the three fields, but we will be interested in the case D — A only 
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and 01 = = Ou and O3 = O^, yielding di^di^d^^ d. Then, Eq. (13) gives 

Ts{puP2; d) = C,J dSi dS, ^^,^^,^,^^,^^,^,^^,^,^, e^^^-^e^^--^ (14) 

for the Fourier transform of the three-point function, where is a constant. Insertion of 
the resolution of the identity in the form 

1 = J d'z5^'\z - {x^ - X2)] - j d\e-^''-[^-(---^)], (15) 

yields 

/i4 
^r2(g; d/2)V2{pi - q; d/2)r2{p2 + q; d/2), (16) 

up to an overall d-dependent constant, which is absorbed into C^. Here 
r2(A;; g) oc (— A;^ — ie)^~^ is the propagator of the scalar field with scaling dimension g as 
in Eq. (6). Since this resembles a standard loop integral, one can introduce Feynman pa- 
rameters and end up with the result 



r3(Pi,P2;rf) = -ie-'^'^C, (^-J Ti(pi/s,p'Js), (17) 
where 

Ti{pI/s,pI/s) = r(2-rf/2)3(4!r)2 / "^^^ "^^^ '^^^^'^^^^ + ^2 + ^3 - 1) 

x(^^J {ymysr'^\ (18) 

-5 = (pi +^2)^ and A = X1X2P2/ s + XiX^pl/ s + X2X3. Now, integration with respect to one of 
the Feynman parameters, say y^, can be performed trivially due to the S function, leaving 
the integrations 

1 1-2/1 
dyi / d|/2 • ■ ■ . 



Next, we perform the following change of variables 

yi = 1 - y2^(^P (19) 

with the Jacobian being equal to a;. The resulting integral is over the unit square and the 
integration with respect to u can be performed analytically, yielding the final result 



2 ' 2 ' ' Ap + B-Bp 



(20) 
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with 2F1 being the hypergeometric function. The integrand is singular at both endpoints, 
but with an appropriate integration method, it can be integrated numerically. The result 
for three different values of d, i.e., 1.1, 1.5, and 1.9, are shown in Fig. 1, which were also 
calculated in Ref. [2] for d = 1.1 and d = 1.9. Just as in the case of the propagator, this is 
the exact three-point function (before coupling to the SM). 

If one requires scale invariance only and not full conformal invariance of the hidden 
sector, the three-point function is not unique. Using self-couplings of continuous mass 
fields, it is possible to obtain nonunique three-point functions with simple expressions in 
momentum space [24] . However, because of the nonuniqueness of these three-point functions, 
this approach is less suitable for phenomenological studies. 

Similarly, the four-point function has the general form [25, 26] 

{o,{x,mx,mx,mx.)) = p ^) n(4.)^<^./6-K+<^.)/2, (21) 

for some function p that cannot be determined by conformal covariance alone. For a single 
field of scaling dimension d, one finds 



P 



\'^13-^24 ■^14-^23/ 



{0{xMx,)0{xs)0\x,)) = J- (22) 



for p' related to p, which, in principle, can be used to study six-body final states. For this, 
one would have to assume a specific form of p' and then perform a Fourier transform of the 
result. This analysis, however, is beyond the scope of the present work. 



IV. PHENOMENOLOGY 



In this section, the main results of this work will be presented. These results will be based 
on an unparticle-induced effect qualitatively different from those of unparticle production 
and effects mediated through the propagator. This is the effect of unparticle self-interactions, 
first introduced by Feng, Rajaraman, and Tu [2], which for example can mediate processes 
with four-body final states at the LHC. While the authors of that paper considered the 47 
final state only and the experimental bound on the unparticle self-interaction strength using 
Tevatron data in the same channel, we will in this work consider also the 2'~f2i and Ai final 
states at the LHC, including those final states in which two of the leptons are neutrinos. In 
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FIG. 1: The function Ti{pl/s,P2/s) for d = 1.1 (top panel), d = 1.5 (middle panel), and d = 1.9 
(bottom panel). The graphs are shown as a function ofpf/s for the following six discrete values: 
pI/s = 10"^ (thick, dash-dotted), 10""^ (thick, dashed), 10"^ (thick, solid), 10"^ (thin, dash- 
dotted), 10~^ (thin, dashed), and 1 (thin, solid). 
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addition, the M channel at the Tevatron will be used to give more stringent bounds on the 
unparticle self-interaction strength and thus also on the LHC cross sections. 

A. Unparticle Self-Interactions at the LHC 

As discussed in Sec. Ill, three- and higher-point functions are in general nonvanishing in 
conformal field theories. These higher-point functions can mediate processes such as 

where Q is any SM particle, Q is its antiparticle, and the number of unparticles are two 
or more. If the conformal sector is strongly coupled, the creation of additional high-energy 
unparticles does not suppress the rate. This is in marked contrast to all SM processes where 
processes with additional high-energy particles are strongly suppressed. 

The natural choice when dealing with these kinds of self-interactions is to concentrate on 
the effects of the three-point function. First, as discussed in Sec. Ill, the three-point function 
is completely specified up to an overall constant by conformal invariance, while the four- 
and higher-point functions are not. Second, this is the lowest order at which the resulting 
signal will have a strongly suppressed background. The unparticle three-point function can 
mediate processes such as 

pp^U ^UU ^ 47, 272^, 

and many others at the LHC. The basic parton-level Feynman diagrams are similar to the 
one shown in Fig. 2. The initial state can be either a ghion or a quark- antiquark pair, while 
the final state consists of two SM particle- antiparticle pairs, i.e., two of the pairs gg, qq, 
W'^W~, ZZ, 77, vV, and HH. In principle, one could also couple SM particles to one 

of the legs in the unparticle three-point function only, while amputating the other leg. The 
amputated three-point function has the form of the ordinary one, but with the amputated 
field replaced with its conformal partner with scaling dimension 4 — 0? [27] . This would lead 
to processes like 

pp^U ^2-fU,2eU 

with a phase space for two particles and one unparticle. Note that, although the addition 
of unparticle lines from higher-point functions leads to no suppression, coupling of these 
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to SM particles does lead to suppression simply because the particle-unparticle coupling is 
assumed to be small. 

In producing the results of this work, the program CompHEP [28, 29] has been used. It 
includes the Monte Carlo integration algorithm VEGAS [30], which is used to integrate over 
phase space, and the CTEQ6L1 parton distribution functions (PDFs) [31]. No high-energy 
physics program that we are aware of has the ability to incorporate such a complicated 
momentum dependence as the one entering through the unparticle three-point function. 
Therefore, in order to use the user-friendly environment of CompHEP when performing 
phase space and PDF integrations, we have altered the C code to be able to enter the values 
of the three-point function. This makes it difficult to perform more complete simulations 
including the interference with SM diagrams and diagrams with one or two unparticle prop- 
agators. However, including the interference with SM diagrams will not be necessary at the 
Tevatron, since they are negligible (for the 47 channel) or can be made so by imposing appro- 
priate cuts (for the 4£ channel) . This also means that the backgrounds will be unimportant 
at the LHC since the SM-unparticle couphngs become more relevant at higher energies. If 
A4 is larger than a few TeV or d is not very close to unity, inclusion of the diagrams with one 
or two unparticle propagators will also not be necessary. The reason is that a diagram with a 
single unparticle propagator contains two SM-unparticle couplings, and thus are suppressed 
by 1/A4'^ in the cross section. In the same manner, two unparticle propagators means a 
I/A4'' suppression. Of course, the contribution from the unparticle three-point function will 
be suppressed by l/A^'', but this can be counteracted by potentially large values of Cd- 

For leptons and photons in the final state, there are six different squared ampfitudes, 
corresponding to the different subprocesses 

99,<N^ 7777, 77^A ^1^1 4^2- 

When calculating observables, the convention is to set Cg — — 1 and e'^{c{Y — 2tt 
[11, 32]. In addition, following Ref. [2], A4 = 1 TeV is chosen. This means that the only free 
parameters we have left are d and Cd- 

It is worth emphasizing that there is no theoretical knowledge on what the values of Cd 
are in typical models of the unparticle sector. Also, there are no model independent bounds 
on Cd coming from more general considerations, although the authors of Ref. [33] beheve 
their method could be used to derive such bounds. As will be shown, generation of the 
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relevant four-body final states requires » 1. Because of the lack of knowledge on typical 
values of C^, one cannot say if this makes it impossible to probe some models of the hidden 
sector through their self-interactions. In addition, if the relevant final states were to be 
observed at the LHC, this could not be used to extract any further information about the 
nature of hidden sector except the scaling dimension itself. 

It is easy to observe that all cross sections will be proportional to C\l and that the 
values of and A4 have no effect other than an overall factor in the cross sections. For 
numerical calculations, the value = 1 is used and can then be given the appropriate value 
to yield the actual cross sections. Thus, the value of A4 chosen has no effect on the total 
cross sections in the above processes if a change in A4 can be accompanied by a change in 
the unparticle self-coupling constant C^- However, this also means that an overall change 
of scale of the dimensionless couplings between the unparticle sector and the SM can be 
compensated by a change in C^. For the sensitivity reach of the LHC on the unparticle 
couphngs to photons and leptons, see e.g. Ref . [34] . Because of the unitarity bound d>\ 
and the more practical bound d <2, only this range on the scaling dimension is considered. 
More specifically, the values d = 1.1, 1.5, and 1.9 are used as representative values, hopefully 
giving a picture of the differences caused by the different scaling dimensions in the chosen 
range. 




FIG. 2: The Feynman diagram for the subprocesses gg/qq ^ W — > UU /i/i/2/2- 
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B. Four-Photon Final States 



The four-photon final state mediated by unparticles has previously been studied in 
Ref . [2] . We include a similar analysis for completeness and also give distributions of invari- 
ant masses, which can be used to identify unparticles as the source as well as identify the 
scaling dimension. Using the Feynman rules given in Sec. II A, one finds that the squared 
amplitudes for the gg — )■ 47 and qq 47 subprocesses are 

2 4 

]m^^ = 2'°^ba ■ PbfiPi ■ P2)\P3 ■ P^f |r3(g2, ^3; d) f , (23) 
2^c^e^(c{)^ 

respectively. These expressions do not include the factor of 1/4! needed to compensate for 

the overcounting caused by the integration over phase space for four identical particles. By 
dimensional analysis, the total partonic cross sections are given by 

\3d 

\ 3d 2 1 

where /| and are dimensionless numbers, which can be extracted from the Monte Carlo 
integration and are listed in Table I for different values of d. 
The final state photons are required to satisfy the conditions^ 

p^>15GeV, \7j'\<2.5 (27) 

with pfp and r/* being the transverse momentum and the rapidity of the i-th photon, respec- 
tively, and i = 1,2,3,4. The resulting reference cross sections o"^^^^'^^^ (i.e., using Cd = 1) 
at the LHC with ^/s = 14 TeV are presented in Table I. Note that our cross sections differ 
slightly from the results presented in Ref. [2], which is most likely due to different parton 
distribution functions. As argued in Sec. IV E, Tevatron measurements give constraints on 
Cd- Multiplying the sum of the reference cross sections by this number squared, yields up- 
per bounds on the LHC cross sections, 174^^'™'^'', which can also be found in Table I. The 

^ All photons have transverse momenta well above 15 GeV. Thus, imposing slightly harder p^-cuts will 
have no effect on the signal. 
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allowed cross sections are observed to depend strongly on d, from only 6.0 fb for d = 1.1 to 
1.2 • 10^ fb for d = 1.9. 

We have also generated a large number of events, enabling the calculation of differ- 
ent distributions. In Fig. 3, the distribution of the two-body invariant photon masses 
rriij = ^ (pi + pj)^ are shown. Since there are six different ways of combining four iden- 
tical particles, each event will contribute six counts. In Fig. 4, the distribution of the full 
four-body invariant mass 1714^ — v^(p7+^2~+^3-i-P4)^ is shown. These spectra are pre- 
sented for the gg initial state only, which dominate at the LHC. However, after having 
generated spectra for a sample of qq initial states and different values of d, one finds almost 
perfect agreement with the spectra for the gg initial state. Thus, if for some reason the 
unparticle field does not couple to gluons at all, the resulting spectra would be virtually 
identical although the total cross sections would be different. What is evident is, first of all, 
the absence of any resonance, owing to the scale-invariant nature of the unparticle sector. 
Second, these distributions depend very strongly on the unparticle scaling dimension. This 
is a general property of all the spectra presented in this work; they all become very much 
softer when d decreases from 1.9 to 1.1. 

The SM background consists mainly of direct diphoton production with two photons 
attached to the quark lines. With the same cuts as above, we find the SM background to 
be approximately 0.1 fb. Obviously, if the unparticle cross sections are anywhere near the 
values allowed by Tevatron data, there is no need to consider the interference effects with 
the SM. Thus, if there was to be an excess of 47 events at the LHC, the given distributions 
could possibly be used to identify unparticles as the source, as well as the correct value of 
the scaling dimension d. Since there is no distinction between the couplings to photons and 
gluons, the distributions of gluons in the Ag and 2g2^ channels would be identical. Although 
the total cross section will be larger in this case due to the multiple colors of the gluons, the 
large QCD background will make the observation in this channel more difficult. 
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FIG. 3: The distributions of photon two-body invariant masses in unparticle-mediated 47 events 
at the LHC for three different values of the scahng dimension d = 1.9 (red sohd curve), d = 1.5 
(black dashed curve), and d = 1.1 (blue dashed-dotted curve). The bin width used is 25 GeV. Note 
that each event contributes six counts. 




'".JGeV] 



FIG. 4: The distributions of photon four-body invariant masses in unparticle-mediated 47 events 
at the LHC for three different values of the scaling dimension d = 1.9 (red solid curve), d = 1.5 
(black dashed curve), and d = 1.1 (blue dashed-dotted curve). The bin width used is 75 GeV. 
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29 


Combined 47 and M bound on Ca 


450 


1.8 • 10"^ 


1.9 • 10^ 


LHC.rcf rr, 1 


2.9 • 10""' 


1.1 • 10-5 


3.2 • 10-5 


LHC,rcf rr, 1 


1.0- 10^6 


3.8 • lO-'^ 


1.1 • 10-6 




6.0 


3.7 • 10^ 


1.2 • 106 


LHC,ref r^, i 


3.8 • 10-^ 


1.1 • 10-5 


1.0 • 10-5 


LHC,ref rr, i 

^55^272^ [ibj 


2.0- 10~^ 


4.0 • 10-^ 


3.7-10-'^ 


LHC.max rn i 

^272/ [fb] 


81 


3.7 • 10^ 


3.7-105 


LHC.ref r^n^l 
<^g5->2e2M [f^l 


2.6 • 10-3 


4.6 • 10-6 


7.7-10-^ 


LHC.ref rrn i 
^qg^2e2M [fb] 


2.7-10-^ 


2.2 • 10-"^ 


2.6 • 10-^ 


LHC.max rri i 
^^2e2/. [fb] 


580 


1.6 • 10^ 


2.9 • 10^ 



TABLE I: The dimensionless constants /J'''^, /i^'"^, and j^'"^, appearing in the partonic cross sections, 
as well as the unparticle 47, 272^, and 2e2/x reference cross sections at the LHC for d = 1.1, 1.5, 
and 1.9. The reference cross sections are calculated with = 1, while the actual cross sections 
are proportional to Cj, which is bounded by Tevatron data (see Sec. IV E), giving upper bounds 
on total cross sections for all final states at the LHC. 
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C. Two-Photon and Two-Lepton Final States 



Simply replacing two of the photons with one charged lepton-antilepton pair, yields the 
subprocesses gg — > 272^ and qq — > 272^ with squared amplitudes 

2 2 2/ f\2 

\M,,-,2,2e\' = (^^ "'"'l^d APa ■ Pi>r{pi ■ P2)^(P3 • P4) iTsiq^, qs', d)\' , (28) 
320^6*^(04)"' 

\Mqq-^2j2ef ^ y ^ ^g/ ^^(Pa " Pb){pi ' ^2)^(^8 ' Pi) \^3{q2, qs', d)f , (29) 

respectively. These expressions should be multiphed with 1/2! to account for the identical 
photons in the final state. Just as in the four-photon case, the partonic cross sections can 
be written as 

/ s \ ^'^ v'^ 1 

^{s)gg-,2'r2e^KCly-j^J JJJ^^^^ (3°) 

\ 3(i / 2 \ 2 -1 

W (t) iTG^f^- P^' 

where the proportionality factors are again listed in Table 1. The same procedure as in the 
four-photon case also yields the reference and maximum cross sections at the LHC, where 
the same single-particle cuts have been imposed. The allowed cross sections are of the same 
order of magnitude as the 47 channel, although the difference between different values of 
the scaling dimension is now smaller. Still, there is the possibility to have very large cross 
sections for large d. 

The Pt distributions of the photon with the highest transverse momentum and the distri- 
bution of the charged lepton invariant masses arc shown in Figs. 5 and 6, respectively. Again, 
one observes the dramatic softening of the distributions as the scaling dimension decreases 
and just as in the four-photon case, almost perfect agreement between the distributions from 
the gg and qq initial states are found. 

Using CompHEP, the SM cross section is calculated to be of the order of 30 fb per lepton 
flavor. There will be resonances at charged lepton invariant masses and Mz — 90 GeV. 
Thus, we demand that the charged lepton invariant masses exceed 110 GcV, which yields 
a reduction of the cross section to approximately 1 fb. As can be observed from Fig. 6, 
this cut will have a very small effect on the signal, especially for d not very close to unity. 
Accordingly, the SM backgrounds will be negligible for anywhere close to its allowed 
values. In analogy with the four-photon case, since the unparticle field is assumed to couple 
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FIG. 5: The pT distributions of the highest-pr photon in unparticle-mediated 272^ events at the 
LHC for three different values of the scahng dimension d = 1.9 (red sohd curve), d = 1.5 (black 
dashed curve), and d = 1.1 (blue dashed-dotted curve). The bin width used is 25 GeV. 
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FIG. 6: The distributions of the charged lepton invariant masses in unparticle-mediated 272^ events 
at the LHC for three different values of the scaling dimension d = 1.9 (red solid curve), d = 1.5 
(black dashed curve), and d = 1.1 (blue dashed-dotted curve). The bin width used is 25 GeV. 



to quarks and gluons, one should also find qqgg, qq'j'j, and 
spectra, but larger cross sections. 



final states with identical 
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D. Four-Charged Lepton Final States 



For four charged leptons in the final state, the squared amphtudes are given by 




■V^iPa ■ Phf{Pl • P2)(P3 • Pi) \^z{(l2, d) 



(32) 





v^{Pa ■ Pb){pi ■ P2){P3 ■ pd \^3{q2, ga; d)f , 



(33) 



3 Af 



respectively. The relevant final states are 4e, 4/i and 2e2/i. The cross sections can again be 
found in Table I and are given for the 2e2// final state. Similar to the previous final states, 
the partonic cross sections are 



In Fig. 7, the pr distributions of the highest-p^ lepton are shown, while Fig. 8 depicts 
the distribution of the invariant masses of three leptons of which two have the same charge. 
Since there are two ways of combining three leptons, each event will contribute two counts. 
Again, there is a strong dependence on d. 

The four-charged lepton channel is a very clean and important channel when it comes 
to the discovery potential of the Higgs boson (H — )■ ZZ — > 4i) at the LHC Thus, one can 
compare the cross sections originating from the unparticle self-interactions to the already 
existing calculated backgrounds in connection with the Higgs boson [35, 36]. With the cuts 
Pt > 5 GeV and \r)\ < 2.5, the irreducible background is about 70 fb, coming mainly from 
two Z-bosons decaying into four charged leptons. Furthermore, with cuts on the invariant 
masses of the charged leptons, it should be easy to reduce this part of the background to a few 
fb. The reducible backgrounds stem from the ti — )■ W^bW^b —fM + X and Zbb M + X 
processes with a cross section of 7.0 ■ 10'^ fb. Since the events coming from these processes 
contain b jets in the final state, they can be reduced below the irreducible levels, using cuts 
on the track isolation and impact parameters without significant reduction of the signal. 
Thus, the possible cross sections from unparticles summed over all channels, 10^ fb — 10^ fb, 
are potentially huge compared to the backgrounds. 




(34) 



(35) 
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FIG. 7: The pT-distributions of the highest-pT lepton in unparticle-mediated 4i events at the LHC 
for three different values of the scaUng dimension d = 1.9 (red sohd curve), d = 1.5 (black dashed 
curve), and d = 1.1 (blue dashed-dotted curve). The bin width used is 25 GeV. 




FIG. 8: The distributions of the three-body invariant lepton masses in unparticle-mediated 41 
events at the LHC for three different values of the scaling dimension d = 1.9 (red solid curve), 
d = 1.5 (black dashed curve), and d = 1.1 (blue dashed-dotted curve). The bin width used is 
25 GeV. Note that each event contributes two counts. 
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E. Tevatron Bounds 

In order to give bounds on the cross sections of processes mediated by unparticle self- 
interactions, one needs bounds on the constant appearing in the three-point function in 
Eq. (17), which determines the strength of the self-interaction. However, there is no known 
way to bound Cd theoretically, and thus, one needs to turn to experiments. These bounds 
could come from many different experiments. Since unparticle effects arc most easily probed 
at high-energy colliders and there are no LHC data available yet, it is natural to turn to 
the existing Tevatron data. It is important to remember that these bounds, and hence, the 
bounds on the LHC cross sections, depend on the values given to the dimensionless coupling 
constants between the SM and unparticle fields. For example, if the unparticle field does 
not couple to gluons at all, the allowed cross section at the LHC will be smaller, since the 
gg initial state is dominating at the LHC, while the qq initial state is dominating at the 
Tevatron. 



1. The 47 channel 

In this section, we repeat the analysis performed in Ref. [2], but find less stringent con- 
straints on Crf. This is because we find smaller cross sections^ and we take into account the 
single photon detection efficiency. 

The SM background at the Tevatron is calculated to be approximately 0.01 fb. Using 
0.83 fb~^ of data, the DO Collaboration found no candidate events [37]. Thus, the 95 % 
C.L. bound equal to three on the mean of a Poisson process [38] becomes the bound on Cd 

' 0.83fb-VJ7'^^fe4' 

where aj^^'^'^^ is the total reference cross section at the Tevatron and e-y = 0.9 is taken as 
the single photon detection efficiency. The values of aj^"^'^^^ and the bounds on Cd can be 
found in Table II. These bounds are approximately 4 times less restrictive in terms of cross 
sections than the bounds found in Ref. [2] . However, this does not matter, since the bounds 
from the 4£ channel calculated in Sec. IV E 2 will be even more restrictive. 



* However, the cross sections are almost identical to the ones in Ref. [2] when we remove all our cuts. 
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d=l.l 


d = 1.5 


d = 1.9 


Tev.ref rn-i 
<^47 [fb] 


8.5 • 10-9 


1.0- 10-1° 


6.1 • 10-12 


Bounds on Cd 


2.6 • 10^ 


2.3 • 10^ 


9.5 • 10^ 



TABLE II: The Tevatron reference cross sections in the 47 channel and the resulting upper bounds 
on Cd- 

2. The 4£ channel 

The DO Collaboration has also searched for four-charged lepton final states [39-41] with 
an approximate integrated luminosity of 1 fb"! in all channels. Since electrons and muons 
arc detected in different ways in the detector, one needs to separate their analyses. 

For the 4/i channel, the relevant geometric cut is \ri\ < 2 and the kinematical ones 
are pt > 15 GeV and cosctj < 0.96, where the Oij's are the angles between all six muon 
pairs. Also, the two invariant masses of the muon pairs in at least one of three possible 
pairings must both be larger than 30 GeV. In addition, there are losses in acceptances 
due to the requirement on having a track match, the quality of the signal of the detected 
muons, and the efficiencies of the triggers. We implement the cuts on rj, pt, cosctj and 
also the more conservative demand that the invariant masses of M12 = a/ {pi + ^2)^ and 
M34 = V (p3 + P4)2 are both larger than 30 GeV. To eliminate the SM background, which 
is dominated by the intermediate states with two Z bosons or photons, the additional cut 
Mi2,M34 ^ [Mz -4rz,Mz + 4rz] = [81.2 GeV, 101.2 GeV] is imposed. We approximate 
the acceptances for the track match, detector quality cuts, and trigger efficiencies to be the 
same as for the ZZ intermediate state in Refs. [39-41], giving a total acceptance loss of 
approximately 64^ = 0.5. 

In the 4e channel, the same requirements on pt and the invariant masses are applied. 
However, the rapidity cuts for electrons are instead I77I < 1.1 or 1.5 < I77I < 3.2. The spatial 
separation R — v^(?77^-^7j)2~-i-"(0i"^-^j)^ is required to be larger than 0.4. In the same way 
as for the A/i channel, the additional acceptance loss is approximated to be 0.5 from the cut 
on R, the electron quality cuts, and the efficiencies of the triggers. 

Finally, for the 2e2/i channel, the same single-particle cuts as in the 4e and 4fi channels 
are applied. Here cos a < 0.96 for the two muons and R > 0.2 for all the final state particles. 
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The same procedure for the 4e channel yields an additional acceptance loss of 0.5. 

The SM background is now approximately zero and no events with our requirements on 
the invariant masses were found. Thus, we obtain the 95 % C.L. bound on 

(J2 ^ 5 

1 fU-l f Tev,ref , Tev,ref , Tev,reA ' 

1 fb e4e (a^^ + a^^ ' + ^2^2^ j 

The cross section in each channel and the calculated upper bounds on Cd can be found 
in Table III. One observes that the bounds from the Ai channels are about 3.3 • 10^, 160, 
and 25 times more severe than the ones from the 47 channel in terms of cross sections for 
d — 1.1, 1.5, and 1.9, respectively. Thus, the 4£ final states produce more restrictive upper 
bounds on the cross sections than the final states. In conclusion, this means that the 
upper bounds on the cross sections from the Ai final states should replace the ones obtained 
earher in the hterature [2]. 





d= 1.1 


d = 1.5 


d= 1.9 


Tev.ref rrui 


4.6 • 10-6 


4.7 • 10-9 


4.4 • 10-11 


Tev,ref rrul 


3.9 • 10-6 


1.8 • 10-9 


1.8-10-11 


Tev.ref rri-i 


2.1 • 10"^ 


1.2- 10"^ 


1.1-10-10 


Bounds on Cd 


450 


1.8 • 10^ 


1.9 - 10^ 



TABLE III: The Tevatron reference cross sections in the M channels and the calculated upper 
bounds on Cd- 



F. Missing Transverse Momentum 

The couplings of SM neutrinos to a scalar unparticle field are usually taken as [42-45] 

_ 

-i-dZiT^aT^nOu, (36) 

where a, ^ = e, ij,,t are flavor indices. For flavor conserving couplings, the matrix A = (A"^) 
is, of course, diagonal. This can lead to the decay of heavy neutrinos to lighter ones [42] 
and other interesting phenomenological consequences. Here, however, the focus will be on 
processes such as 

pp — >■ 27i/I7, 2^i/I7. 
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The couplings in Eq. (36) will lead to the same squared amplitudes as given in Eqs. (28) 
and (32), up to an overall factor. 

One cannot say anything about the total cross sections for these kinds of processes unless 
one has knowledge about the couplings A°'^. For A"'^ of the order of unity and A3 of the 
order of 1 TeV, the resulting cross sections would be of the same order of magnitude as 
the ones calculated for the 272^ and 4£ channels in Sees. IV C and IV D, respectively. The 
distributions of missing transverse momentum, for 2'y2v and 212^ events at the LHC are 
shown in Figs. 9 and 10, respectively. Compared with the events expected from SM processes 
[46, 47], the events mediated by unparticles are observed to contain large f^, especially for 
larger values of d. In addition, since the total rates could be relatively large, a high cut on 

could be used to search for unparticle physics in the 27 +frp and 21 + '^^ channels as well. 

V. SUMMARY AND CONCLUSIONS 

Unparticle physics is a relatively new idea for an extension of the SM of particle physics, 
introduced by Georgi in 2007 [1] . The basic idea is that there might exist a hidden conformal 
sector, which is coupled to the SM through effective operators. Thus, unparticle physics is 
the study of the phenomenological consequences of the SM interactions with such a hidden 
sector. 

In this work, the implications of the unparticle self-interactions, entering through the 
three-point correlation function, for LHC phenomenology have been investigated. For weak 
enough couphngs of the unparticle sector to the SM, i.e., for a high enough unparticle 
scale A4, unparticle physics will be impossible to observe in all previously studied processes, 
i.e., those including unparticle propagators and external lines. Because of the new degree 
of freedom entering through the self-interactions, this does not necessarily apply to the 
processes studied in this work, assuming large values of are allowed. Thus, it is not only 
possible that these processes could be used as a favorable way to look for unparticle physics, 
but also that it might be the only way. 

The final states examined were the 47, 272^, and M final states, the first of which has 
been studied previously in the literature [2]. In addition, we have used the analysis of 
the 4£ channel at the Tevatron to calculate constraints on the constant appearing in the 
unparticle three-point correlation function, enabling determination of the maximum allowed 
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FIG. 9: The missing transverse momentum in unparticle-mediated 271^!/ events at the LHC for 
three different values of the scahng dimension d = 1.9 (red solid curve), d = 1.5 (black dashed 
curve), and d = 1.1 (blue dashed-dotted curve). The bin width used is 25 GeV. 
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FIG. 10: The missing transverse momentum in unparticle-mediated 2ii/T' events at the LHC or 
three different values of the scaling dimension d = 1.9 (red solid curve), d = 1.5 (black dashed 
curve), and d = 1.1 (blue dashed-dotted curve). The bin width used is 25 GeV. Note the different 
scale on the pj, axis as compared to Fig. 9. 

cross sections at the LHC. The results take the form of upper bounds on the LHC cross 
sections and distributions of transverse momenta, invariant masses, and missing transverse 
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momenta. 

The allowed cross sections have been found to depend strongly on the scaling dimension d 
of the unparticle sector. In general, they are much smaller for d close to 1 and increase as d 
increases towards 2 in which case the cross section can be as large as 10^ fb. For d = 1.9, we 
obtained 10® fb for the 47 final state, 10^ fb for each of the the 272^ final states, and 10^ fb 
for each of the M final states. However, for smaller values of 0?, the allowed cross sections 
are much smaller due to our new constraints from the M channel at the Tevatron, which 
are indeed much stronger and replace the upper bounds from the 47 channel. In addition, 
the calculated distributions have been observed to depend strongly on the scaling dimension 
and become much softer when it decreases towards 1. Finally, if there was to be an excess of 
events of the relevant final states at the LHC (possibly enormous), unparticle physics would 
be one possible description. The computed distributions could then potentially be used to 
identify unparticles as the source, as well as the correct value of the scaling dimension d. 
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